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HYPERBOLIC CONSERVATION LAWS 
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Abstract. This paper is concerned with the initial-boundary value problem for a nonlinear hyper- 
bolic system of conservation laws. We study the boundary layers that may arise in approximations 
of entropy discontinuous solutions. We consider both the vanishing viscosity method and finite 
difference schemes (Lax-Friedrichs type schemes, Godunov scheme). We demonstrate that differ- 
ent regularization methods generate different boundary layers. Hence, the boundary condition 
can be formulated only if an approximation scheme is selected first. Assuming solely uniform L°° 
bounds on the approximate solutions and so dealing with L°° solutions, we derive several entropy 
inequalities satisfied by the boundary layer in each case under consideration. A Young measure is 
introduced to describe the boundary trace. When a uniform bound on the total variation is avail- 
able, the boundary Young measure reduces to a Dirac mass. Form the above analysis, we deduce 
several formulations for the boundary condition which apply whether the boundary is character- 
istic or not. Each formulation is based a set of admissible boundary values, following Dubois and 
LeFloch's terminology in "Boundary conditions for nonlinear hyperbolic systems of conservation 
laws", J. Diff. Equa. 71 (1988), 93-122. The local structure of those sets and the well-posedness 
of the corresponding initial-boundary value problem are investigated. The results are illustrated 
with convex and nonconvex conservation laws and examples from continuum mechanics. 



1. Introduction. 

This paper considers the initial-boundary value problem for an hyperbolic system of con- 
servation laws 
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supplemented with 

(1) an initial condition at time t = 



u(x, 0) = ui(x), 



x > 0, 



(1.2) 



(2) 




(1.3) 



(3) 



u(0,t) = u B (t), 



t > 0. 



(1.4) 



Indeed the hyperbolic problem (1.1)-(1.4) is usually not well-posed when the boundary 
data is required to be assumed in the (strong) sense (1.4), even when (1.1) is a linear system 
(cf. Kreiss [28]). It is the objective of this paper to provide a general framework which leads to 
(mathematically correct) formulations for the boundary condition. Following Dubois-LeFloch 
[15], our strategy is to reformulate (1.4) in the (weak) form 



where S(uB(t)) C U is a time-dependent set (the set of admissible boundary values) to be 
defined from the boundary data, and u(0+,t) is the trace (its existence is discussed in this 
paper) of the solution u at the boundary. We shall consider several methods of approximation 
for the problem (1.1)-(1.4), including the artificial vanishing viscosity method and a class of 
finite difference schemes, for which the boundary condition (1.4) can be easily implemented. As 
the approximation parameter goes to zero, a sharp transition layer generally develops near the 
boundary {x = 0} and the limiting solution does not satisfy the boundary condition (1.4). Our 
aim in this paper is to provide some contribution to the following program: perform a rigorous 
analysis of the boundary layer for weak solutions, then derive several suitable definitions for the 
set in (1.5), and finally investigate the structure of the latter to decide whether the boundary- 
value problem is well-posed. 

In (1.1), U is assumed to be a convex and open subset of R N , the flux- function / : U — > R N 
to be a smooth mapping, and the initial data uj to belong to L°°(R + ,U). It will be convenient 
to assume that the boundary data ub has bounded total variation on any interval [0, T] for all 
T > 0. It is assumed that (1.1) admits at least one strictly convex entropy pair. By definition, 
a pair of functions {U, F) :U — > RxR of class C 2 is called a convex (or strictly convex) entropy 
pair iff VF T = VU T V f and the Hessian matrix V 2 U is non-negative (or positive definite). 
The existence of at least one strictly convex entropy pair implies that (1.1) is hyperbolic. For 
background on hyperbolic systems, we refer to Lax [29, 30, 31], Dafermos [11] and Smoller [45], 
concerning the theory of existence of entropy solutions to the pure Cauchy problem, to Glimm 
[21] and Liu [40] for initial data with small total variation, and DiPerna [12,13] for systems of 
two equations with L°° initial data. 

This paper contributes to establishing a framework for the initial-boundary value problem for 
(1.1). It is intended to pursue the efforts initiated in recent years on this problem (Cf. review 



u{0+,t) E £(u B (t)), 



t > 0, 



(1.5) 
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below). In particular we built upon the recent contributions in Gisclon-Serre [20] and Xin 
[49], who studied the boundary layers associated with the vanishing viscosity approximations 
assuming the solution to the hyperbolic problem be smooth. A formal asymptotic expansion 
is introduced in [20, 49] and the convergence including L 2 error estimates is proven for the 
boundary layer in the smooth regime. For linear hyperbolic systems, Joseph [25] constructed 
boundary layers explicitly and obtained error estimates in L 2 Sobolev space. 

One of the motivations here is to treat several approximation methods simultaneously and 
compare the results obtained with each of them. We consider the vanishing viscosity method, 
a class of Lax-Friedrichs type schemes, and the Godunov scheme. 

In Section 2, we rigorously derive conditions satisfied by the boundary layer, which take the 
form of a family of boundary entropy inequalities and a boundary layer equation. The regularity 
of the relevant traces at the boundary are discussed. The whole analysis is performed by 
assuming only a uniform L°° bound on the approximate solutions; in particular no assumption 
is required on the regularity of the limiting solution to (1.1). Since high frequency oscillations 
in the approximate solutions can not be a priori excluded, the conditions above are formulated 
in terms of a boundary Young measure associated with the boundary layer. Note that, in 
the derivation of Section 2, the boundary is possibly characteristic, i.e. the eigenvalues of the 
matrix Vf(u) may vanish for certain values of u. 

Observe also that, in general, the equations and inequalities we derive depend upon the 
approximation method in use. Fundamentally the boundary condition can not be formulated 
from the mere knowledge of the function ub, but depend upon the underlying "physical" 
regularization. This feature arises in weak solutions to many nonlinear hyperbolic problems. 
See, for instance, the review paper by LeFloch [33] on regularization-sensitive shock waves. 

In Section 3, we introduce several sets of admissible boundary values and investigate their 
local structure. When the boundary is non-characteristic, we establish that the sets based 
on the boundary layer equations are manifold with the "correct" dimension. That is, the 
corresponding initial-boundary value problem is well-posed, at least for constant boundary 
and initial data (a generalization to the Riemann problem). We also prove a similar (but 
stronger) result for the set based on the boundary layer equation derived by the Godunov 
scheme. Strictly speaking this scheme does not produce any boundary layer; however analyzing 
that scheme leads to a formulation of the boundary condition as it was first pointed out in 
[15, 16]. We recall that setting the boundary condition via an upwinding difference scheme is 
a classical idea in the computing literature. 

Sections 4 is ddevoted to studying several examples of particular interest. It is expected 
that, in general, different approximation method for (1.1) leads to a different set in (1.5). 
However we prove in Section 4, for both convex and non-convex conservation laws, that this is 
not the case when N = 1. In other words the boundary layer for the scalar conservation laws is 
independent of the approximation method. The same is true of the linear hyperbolic systems; 
and we conjecture that this also holds for the nonlinear systems in the class with coinciding 
shock and rarefaction curves introduced by Temple [48]. We also consider examples from 
continuum mechanics, i.e. the system of nonlinear elasticity and the system of gas dynamics. 

To complete this presentation, we give a short overview of the literature on the boundary 
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conditions for (1.1). Most of the activity was restricted to scalar equations, i.e. N = 1. The 
pioneering work by Leroux [35] and Bardos-Leroux-Nedelec [4] based on the vanishing viscos- 
ity method provides a derivation of "the" correct formulation of the boundary condition for 
multidimensional scalar conservation laws. Specically, [4] shows that (1.4) should be replaced 
by the weaker statement: 

(sgn(u(O+,*)-A0-sgn(u B (*)-A0) (/(«(0+, t)) - f(k)) >0 for all k E R, (1.6) 

where sgn(a) = — 1 if a < 0, sgn(a) = if a = 0, and sgn(a) = 1 if a > 0. The convergence of 
finite difference schemes, again for scalar equations, is established by Leroux in an unpublished 
work: it is remarkable that the finite difference scheme approach leads to the same formulation 
(1.6) of the boundary condition. The condition is used by LeFloch [32] in order to extend Lax's 
explicit formula [30] to the initial-boundary value problem. Joseph [24] used the vanishing 
viscosity method and the Hopf-Cole transformation to extend Lax's formula for the inviscid 
Burgers equation. Another derivation is given by Joseph and Veerappa Gowda [27]; see also 
Gisclon [18] and LeFloch-Nedelec [34] .We also refer to the paper [47] by Szepessy for a very 
general result of existence and uniqueness. 

The statement (1.6) is a special case (when applied to Kruzkov entropies) of a more general 
inequality: 

F(u(0, t)) - F(u B (t)) - VU(u B (t))(f(u(0, t)) - f(u B (t))) < 0, (1.7) 

which has to hold for every convex entropy pair (U,F). The latter was derived formally 
using the vanishing viscosity method in Dubois-LeFloch [15], who pointed out that (1.7) holds 
even when N > 2 and introduced the notion of set of admissible boundary values, cf. (1.5). 
These inequalities were obtained independently by Bourdel-Delorme-Mazet [8] based on an 
analysis of the characteristics of the system (1.1), and by Benabdallah [5] for a specific system. 
The first result of existence for the initial-boundary value problem for a system was given 
by Benabdallah-Serre [6, 7]: the vanishing viscosity method applied to the p-system of gas 
dynamics converges to a solution to (1.1) satisfying the set of inequalities (1.7). 

The Glimm scheme with various type of boundary conditions was studied by Liu, for in- 
stance [37, 38, 39]. In the case that the boundary is assumed to be non-characteristic and the 
number of boundary conditions is equal to the number of positive eigenvalues of the matrix 
V/, Goodman proves the convergence of the Glimm scheme in his unpublished thesis [22]; 
cf. also Dubroca-Gallice [17] and Sable- Tougeron [42, 43] 

More recently Amadori [1, 2] used the formulation in [15] and proved the convergence of 
a front tracking scheme in the characteristic case. In particular, Amadori establishes that a 
condition of the form (1.5) can be satisfied pointwise except at countably many times. 

Finally we refer to the IMA report [26] by the authors for an extended version of the present 
article. 

2. Boundary Layers in Weak Solutions. 

In this section, we consider sequences of approximate solutions to the initial boundary value 
problem (1.1)-(1.4), and aim at characterizing their limiting behavior near the boundary. 
Here we rigorously derive entropy inequalities satisfied by the boundary layer. We deal with 
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a sequence of L°° functions with uniformly bounded amplitude. As is well-known, for general 
systems of conservation laws, proving the strong convergence of a sequence of approximate 
solutions is an open problem. It seems therefore natural to formulate those entropy inequalities 
in terms of a Young measure (for instance Ball [3] for this concept) associated with the sequence 
of approximate solutions. Further analysis can be performed on a case by case basis only. 

In the following, certain averages will be shown to belong to the space BV(R + ) of functions 
of locally bounded total variation, i.e. measurable and bounded functions w : R + — > R whose 
distributional derivative is a bounded Borel measure on every interval (0, T) for all T > 0. 
We denote by TV^f(w) the total variation, and by HHIw^T) = ||H|l°°(o,t) + TVq (w) the 
norm, of a BV function w on an interval (0,T). By convention, a BV function will be always 
normalized by selecting its right continuous representative. 

2.1 Vanishing Viscosity Method. Let u e be the approximate solutions obtained by solving 



the following parabolic regularization of (1.1)- (1.4): 

d t u c + d x f{u e ) = edl x u% x > 0, t > 0, (2.1) 

u £ (x,0) = u € j(x), x>0, (2.2) 

u e (0,t) = u e B (t), t>0. (2.3) 



The smooth functions uj G L°°(R + ) and u B G BV(R + ) are chosen to be uniformly bounded 
and a.e. convergent approximations of the corresponding data ui and ub- We assume the 
existence of a (smooth enough) solution u e to the problem (2.1)-(2.3). Note that compatibility 
conditions at (x,t) = (0,0), such as wf(0) = u B (0), are implicitly required. We shall also 
assume that 

u e is uniformly bounded in L°°(R x R+). (2.4) 
We introduce a new function v e by setting 

v e (y,t) = u e (ey,t), (2.5) 

so that the system of equations (2.1) transforms into 

ed t v e + d y f(v e ) =d 2 yy v\ (2.6) 

It is expected that the (e — > 0) limit of the v e 's will give us a good description of the boundary 
layer at x = 0, at least under additional assumptions, although a different scaling may more 
adapted in certain circumstances. 

By definition (e.g. Ball [3]), a Young measure associated with a sequence u e satisfying 
(2.4) is a weak-star measurable mapping v from the (x,t) plane to the space Prob(W) of all 
probability measures (i.e. non-negative measures with mass one) with the property that for 
every continuous function g :U — > R 

g(u £ ) ^<v,g> weakly-* in L°°(R 2 + ). (2.7) 

In view of (2.4), the functions v e also are uniformly bounded in L°°(R x R+)). We denote by 
/it a Young measure associated with the functions v e . 
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Theorem 2.1. The following statements hold for all convex entropy pairs (U, F) associated 
with the system (1.1), all functions 9 G BV(R+), and any bounded interval (Ti,T2). 
1) When 9{t) > 0, the distribution 

f T2 d f T2 

y^ / < fjL Vtt ,F > 6(t)dt- — < n yjt ,U > 6(t)dt 
Jt x dy J Tl 

is in fact a function of locally bounded variation and thus is defined pointwise as a right con- 
tinuous function. There exists a Young measure /Uo,t, such that the following limit exists and 
is given by /x ,t •' 



pT2 pT2 

lim / < Hy jt ,U > 9{t)dt = <no tU U> 6(t)dt. 

y^0+ J Ti J Ti 

When 9(t) > 0, the function 



/■T2 

x i-> / < u Xit ,F > 6{t)dt 



has locally bounded variation. There exists a Young measure vqj> the "trace" of u Xjt at x = 0, 
such that the following limit exists and is given by u 0jt : 

lim / < u Xtt , F > 6(t) dt = f <v , u F > 6(t)dt. 

x^0+ J Ti J Ti 

When (U, F) = (uj, fj), 1 < j < N, all of the results above still hold when the function 9 has 
no specific sign. 

2) For all < y\ < y2 and in the sense of distributions for t G R+, one has 

F(u B ) + VU(u B )(< v Q ,uf> -f(u B )) ><fi Vl ,t,F> -d y < n yuU U > 

>< Hy a , t ,F> -d y < fiy 2 , t ,U > (2.8) 
>< VQ, t ,F> . 

3) Moreover one has 

t*o,t = $u B {t) a.e. t G R+ (2.9) 

and, when 9 > 0, 

limf/ < n y ,t, F > 9{t) dt - 4- I < l*y,t,U > 9(t)dt] > [ < u , t , F > 9{t) dt. 

y^oo \J Ti dy J Tl J J Tl 

(2.10) 
□ 



K.T. JOSEPH AND P.G. LEFLOCH 



7 



Theorem 2.1 provides a rigorous basis to the formal asymptotic expansion approach. We 
collect here several important remarks, including the property that the Young measures v and 
fi reduce to Dirac masses when a uniform total variation bound is available. 

First of all, the inequalities (2.8) actually hold in the (stronger) sense: 

(F(u B (t)) + VU(u B (t)) (< u , t , f > -f(u B (t))^J 9{t) dt 

T 2 , ,T 2 



> JJ <fji yi , t ,F> e{t)dt-^-{JJ <fjL V , t ,u > e(t)dt)\ y=yi 

> [ 2 <n y2>u F> 9(t)dt-—(f 2 <n v ,t,U > 0(t)dt)\ 

> [ 2 < v ,t, F > 9{t) dt 



V2 



for all non- negative 9 G BV(R+) and all < y\ < x/2- Observe that this is a stronger 
statement than the convergence in the sense of distributions since 9 is a function of bounded 
total variation, not necessarily having compact support in (Ti,T 2 ), rather than a smooth 
function with compact support. All the formulas to be derived in this section hold in this 
sense. Note also that (2.10) is an immediate consequence of (2.8) by taking y — > oo. 
The following inequalities, rigorously derived in Theorem 2.1, 

F{u B ) + VU(u B ) (< v , f > -f(u B )) > < is , F > (2.11) 

will be referred to as the boundary entropy inequalities. They do not refer explicitly to the 
boundary layer itself but only to its limiting values. 

The inequalities (2.8) also contain constraints for the boundary layer. In particular, using 
the trivial entropies (U, F) = ± fj(u)), 1 < j < N, in (2.8) leads us to the equation 

-a y < ^, id >=< ^ ,t, / >, (2.12) 

where the right hand side is independent of the variable y and only depends on t. 

For scalar equations and when the method of compensated compactness due to Murat- Tartar 
applies (i.e., mainly, for systems of two conservation laws), it is known that v is a Dirac mass 
concentrated at a point u(x, t) which is an entropy weak solution. In those two situations, it 
is conceivable that the Young measure \i also would be a Dirac mass. 

If one assumes that fx is a Dirac mass, say 

A*y,f = $v( y ,t) for almost every (y, t) (2.13) 

with v G then the formulas in Theorem 2.1 take a much simpler form. Namely if (2.12) 
holds, then (2.12) becomes what will be referred to as boundary layer equation: 



f(v) -d y v =< i/o,/ > 



(2.14) 



8 



BOUNDARY LAYERS IN WEAK SOLUTIONS 



This is nothing but the equation that would be obtained formally by plugging an asymptotic 
expansion of the form u e (x, t) = u(x, t) +v (x/e, t) + 0(e) in the equations (2.1). More generally, 
if (2.12) holds, the inequalities (2.8) become 

F(u B ) + VU(u B ){< u J > -f(u B )) > F(v( yi ))-d y U(v) ly=Vl 

> F(v(y 2 ))-d y U(v) ly=y2 

> < v ,F > . 

When vq also is a Dirac mass for a.e. t, say z/o,t = o~u (t), f° r instance when u has bounded 
variation in x and so admits a trace at x = in a classical sense, then the boundary layer 
equation (2.14) becomes 

f(v)-d y v=<u J>. (2.15) 
and the boundary entropy inequalities (2.11) take the form 

F(u ) - F(u B ) - VU(u B ) (/(mo) - f(u B )) < 0, (2.16) 

which was derived in Dubois-LeFloch [14, 15] by assuming a uniform BV bound on the u e . 

Note finally that the behavior of (i y j as y — > oo is controled by the set of inequalities (2.10), 
only. If it is assumed that v has a limit in a classical sense and d y v(y,t) — > as y — > oo, then 
we can set 

Uoo(t) = lim v(y,t) 

and (2.10) becomes 

^(foo) > F(uo) for all entropy flux F (2.17) 
(the flux F must be associated with a convex entropy). In fact (2.17) need not imply 

v oo (t)=u (t). (2.17') 

However (2.17) does imply 

f(v oo (t)) = f(u (t)) 

so, in the non-characteristic case i.e. when V/ is invertible, (2.17) implies (2.17'). In the 
characteristic case, (2.17') may very well be violated. This difficulty is related to the choice of 
the scaling in the definition of the functions v e . Cf. the examples in Sections 4 and 5. 

Uniform bounds on the total variation of u £ are available for scalar equations, linear systems 
and systems in the so-called Temple's class having coinciding shock and rarefaction curves. In 
the general case we have : 

Corollary 2.2. Assume that the solutions to the boundary-value problem (2.1)-(2.3) addi- 
tionally satisfy the bound 



/•OO 

TV(u e (t)) := / \d x u e (t)\ dx < C, 
Jo 
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where C is independent of e. Then the Young measures \i and u Q in Theorem 2.1 reduce to 
Dirac masses, i.e., 

fJ> y ,t = Sv( y ,t) and z/ ,t = 5 Uo (t)> 

and the functions u = v(y,t) and -uo = «o(£) satisfy the conditions (2.14) and (2.16) for almost 
every y, t. □ 

Proof of Theorem 2.1. We decompose the proof into several steps. For the whole of this 
proof, we denote by (U, F) a given convex entropy pair. 
Step 1: Preliminaries. 

We gather here several properties of v and fi that are readily obtained. Let us multiply the 
equation (2.6) by the gradient of U and obtain 

ed t U(v e ) + d y (F(v*)-d y U(v*)) = -V 2 U{v*) ■ {d v v',dyV*) 

< 0. 1 ' ' 

Using the definition of the Young measure fi, it is a simple matter to pass to the limit in the 
inequality (2.18). For any 9 E BV and uniformly in y G we have 



/ 2 d t U(v e )9dt\ < | / eU(v e )d t 9dt\ + e\[U(v e ) 9]^ 2 \ 
<O(l)||0|| B v||tf(t; e )IU- - 0, 



e 



so we obtain 



dy^f 2 <Vy,t,F> 9dt-^-J 2 <Hy jt ,U> 9d?j < 0, (2.19) 

which provides the second inequality in (2.8). Therefore time-averages of the function < 
fiy jt , F > —d y < fj,y jt , U > axe noil- incr easing , and so have bounded variation on any compact 
set. The limits as y — > 0+ or y — > +oo exist, although at this stage of the proof, we can not 
exclude that those limits could be ±oo. We shall see later that actually < ^ y ,t,F > —d y < 
fi yj t, U >E L°°. Moreover the function 

/ < fJLy, t , U > 9{t) dt 

has a trace at y = 0, which defines < ^0,4, U >. Note also that (2.19) with the choices 
(U, F) = ±( Uj , fj), 1 < j < N, leads us to 

< Hy jt , f>-8y< flyj, id >= C *(t) , (2.20) 

where C*(t) has to be determined. In fact it will be immediate from the results in Step 5 below 
that 

C*(t) =< ^o,t 5 id > for a.e. t > 0. 
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Similarly, following DiPerna [13] and using the Young measure v x<t associated with u e , one 
can pass to the limit in (2.1) and obtain the entropy inequality: 

8 t < v XtU U>+8 X < v Xtt , F >< 0. (2.21) 

From (2.21), we deduce first that, for any smooth function 9(t) > 0, 

^ jf 2 < v x ,t, F > 9(t) dt< J J < u Xtt , U > d t 9(t) dt < 0(1) \\6\\bv (2-22) 

For 9 fixed, the right hand side of (2.22) is a constant, thus its left hand side is a locally 
bounded Borel measure and the function 



pT2 

g e (x) = / < u Xit ,F> 9(t)dt 



has bounded total variation. Therefore the trace fo,t introduced in Theorem 2.1 exists, at least 
on entropy fluxes. This gives a meaning to the last term in the right hand side of (2.8). In fact 
it is possible to establish the estimate 

TV (g ) < 0(1) \\9\\ BV 

for arbitrary functions 6 E BV. (For such 6, (2.22) can be obtained directly from (2.1).) Thus 
the trace z/ ,t exists for 9 G BV as well. 

Observe that the traces jiio,t and z/Q,t are uniquely determined on entropies and entropy 
fluxes, respectively. They can be easily extended as Young measures defined on the whole set 
of continuous functions, in a non- unique way however. Namely, to construct /Uo,t, take any 
sequence — > and consider a Young measure associated with the sequence of measures 

This completes the proof of the part 1) in Theorem 2.1. 

Step 2: A General Identity. 

It remains to analyze the behavior of \i at the end point y = which shall provide us 
with the desired boundary entropy inequality. We are going to use a general identity which 
immediatly follows from the Green formula applied to (2.6). 

Let 9(t) and <p(y) be smooth functions not necessarily having compact support. We multiply 
the equation (2.6) by VU(v e ) 9(p and integrate over the domain (j/i, y 2 ) x (T 1? T 2 ). Integrating 
by parts and re-ordering the terms, we obtain the identity 

E\ + Ejj + Ej n = Ej v (2.23) 

with 

= -e / U(v e )d t 9ipdydt + e9(T 2 ) U(v e (T 2 )) ipdy - e0(Ti) / U^^)) <pdy, 

(2.24.1) 



Ej = 
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Eh = ~ / F(v e )6d v <pdydt + <p(y2) / (F(v e (y 2 )) - d y U(v e )\ y=y2 )0 

<p{ Vl ) [ 2 (F(v*( yi ))-d y U(v% =yi )9dt, 

U(v e )dd yy <fdydt + d y <p(y 2 ) [ * ' U(v e (y 2 ))9dt 

/ a u(v'( yi ))6dt, 
Jt 



e iii = ~ 



T 2 t-yi 



(2.24.III) 



and 

Ej v = - \ I VU(v e ) ■ (d y v\d y v e )e^dydt. (2.24.IV) 

JT-y J Vl 

In case that 9 > and </? > and since U is assumed to convex, one has 

Ej v < 0, (2.25) 
so we can focus attention on estimating the terms Ej, Ejj and E\ n . 

Step 3: Viscous Flux at the Boundary. 

We prove here that the viscous flux at the boundary, i.e. the function d y v e (0, t), is uniformly 
bounded in a certain sense and we determine its weak limit as e — > 0. We use the identity 
(2.23)-(2.24) with the following choice of parameters: 

supp 6 C [T l5 T 2 ], supp y C [0, 1), = 0, y 2 = 1, (tf, F) = ( Uj , /,), 1 < j < N. 

For <£> fixed, we obtain 

|£fl < O(e)||0|| w , 

r-T 2 /.l r T 2 



Eh = ~ [ 2 f f(v<)ed v <pdydt-<p(0) [ 2 \f(u%)-d y v € (0,.))edt 

JT X JO JTx 

= 0(1) 11^11^-^(0) / 2 (f(u%)-d y v*(0,.))ddt, 



and 



*T 2 /.l /.T 2 



E\ n = -\ \ v e 0d yy (pdydt - u e B 9d y (p(0) dt 
Jti Jo Jti 



IT! JO JT X 

= 0(1) ||0|| L oo. 



Since in this case Ej v = and choosing (p so that ip(0) ^ 0, it follows 

fT 2 

IT! 



I / * (f(u%)-d y v e (0,.))6dt\ < 0(1) + O(e) ||0||b V . (2.26) 

Jt 
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More precisely we can pass to the limit in the identity (2.23) and get 

'Ti 



V?(0) lim / 2 (f(u B )-d y v e (0,t))9dt 

/•T 2 pi pT2 pi PT2 

— / < fx, f > 9d y (p dydt — / / < fx, id > 9d yy cpdydt — d y (p(0) / 

./Ti io JT X JO JT X 



On the other hand, it has been observed in Step 1 that (2.20) holds and < fx, id > has a trace 
at y = 0. Thus one has 

T 2 />1 /-T2 pi 

< fx, f > 0dy(fdydt + / < fx, id > 9d yy (f dydt 



Ti JO JTi JO 

T 2 /•! /.T 2 



= / / C*(t)6d y ipdydt- < fx , id > 9d y <p(0) dt 
Jt x Jo Jt x 

= -/ C*(t) 9<p(0) dt - < fx ,id > 6d y <p(0)dt 

JTl JTl 



and therefore 

"T 2 



p(0) lim / \f(u B ) -d y v e (0,t))9dt 

/•T 2 PT2 PT2 

<p(0) C*(t) 6 dt + d y <p(0) < li , id > 9 dt - d y <p(0) u B 0dt. 

JTl JTl JTl 



Choosing two test-functions (p, one such that (p(0) = but d y (p(0) 7^ 0, and the other such 
that (p(0) 7^ but d y (p(0) = 0, we deduce from the above formula that 



lim ! 2 (f(u B ) -d y v £ (0,t))8dt = j ' C*{t)9dt 

UJr ■'' (2.27) 

/ < fx , id > 9 dt = u B 9dt. 

JT1 Jti 



IT! JT X 

The first statement in (2.27) is the desired convergence result. The second statement is a first 
step toward proving (2.9). 

Step 4: Boundary Entropy Inequalities (I). 

Using (2.27), we are now able to obtain the boundary entropy inequalities. We use the 
identity (2.23)-(2.24) with 

9 > 0, supp 9 C [Ti, T 2 ], <p > 0, supp <p C [0, 00), y x = 0, y 2 > 0, 
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and (U, F) arbitrary. We obtain 

\Ej\ < 0(e) \\0\\bv, 

T 2 ry 2 pT 2 



rU2 

/ F(v € )Od y <pdydt-<p(0) / {F(u%) — d y U(v e )y=o) 9dt 

JO JT X 

I F(v £ )9d y pdydt-p(0) J \F(u B ) - VU{u%) d y v e (0, .)J9dt 



Ti 

pT 2 ry 2 rT 2 / \ 

^-J t J <^F> 6d y <pdydt-<p(0) \F{uB)-VU{u B ){f{uB)-C*{.))j9dt, 
where we have used (2.27) and the fact that u B G BV converges strongly to ub £ BV, and 

fT 2 ry 2 rT 2 

Ej n = - / U(v e )9d yy ipdydt- \ U(u%) 9d y (p(0) dt. 
Jt x Jo Jt x 

Since Ej v < we pass to the limit in (2.23) and get 

„T 2 



p(0) jT 2 (f(u b ) - VU(u B ) (f(u B ) - C*(t)f) 9dt, 

rT 2 ry 2 

-~ / {< V y ,t,F > d y <p+ < /jy, t ,U > dyyip) 9dydt 
Jt x Jo 

pT 2 

+ (fM / (< Vy 2 ,t,F > -8y < Hy,t,U >y = y 2 )9dt 

rT 2 rT 2 

+ d y ip(y 2 ) < Hy a , t ,U > 9dt- d y <p(0) U(u B )9dt. 
Jti Jti 



On one hand, using the test- function <p(y) = 1, we deduce that 

j 2 (F(UB)--VU(u B )(f(UB)-C*(t))j9dt > j 2 (< ^F>+8y< fl,U>y=y 2 )edt 

(2.28) 

which proves the first inequality in (2.8). 

On the other hand, using the function <p(y) = y, we obtain 

fT 2 ry 2 rT 2 

0>- / / < fJLy, t ,F > 9dydt+ y 2 \ (< Hy 2 ,u F > -d y < fi y ,t,U > y=y2 ) 9 dt 
Jt x Jo Jt x 

pT 2 pT 2 

+ / < Vy 2 ,t,U > 9dt- / U(u B )9dt, 
Jt x Jt x 

which as y 2 — > yields 

pT 2 pT 2 

/ U(u B )9dt> lim / </j yt7 U>9dt. (2.29) 
JTi y^°+ Jti 
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In particular, plugging (U, F) = (uj, fj), 1 < j < N, in (2.29), we recover the second statement 
in (2.27), which used together with (2.29) for any fixed, strictly convex entropy U gives: 

\ < no it , U - U{u B ) - VU(u B )(id - u B ) >0dt 

/■T2 

= lim / < Hy t ,U -U(u B ) -VU(u B )(id-u B ) > dt 
< f U(u B )6dt-J U(u B )6dt 



= 0. 



But the function u — > U(u) — U(u B ) — VU(u B )(u — u B ) is positive everywhere except at u B 
where it achieves its global minimum value. It follows that fio }t is a Dirac mass concentrated 
at u B . That proves (2.9). 

Step 5: Boundary Entropy Inequalities (II). 

We now establish the third inequalities in (2.8). We use once more the identity (2.23)-(2.24) 
with now 

6 > 0, supp 9 c [Ti, T 2 ], if > 0, supp cp c [yi, 00), y 1 > 0, y 2 = 00, 
with a function <p depending on e, that is 

<P e (y,t) = <p(ey,t) 
with dp fixed. In that situation one can check that 



/■T2 poo 

Ej = - / U(u £ )d t 9^dxdt 

JT\ Jey 1 

/(•T-2 poo 
Ti io 



-> - / / < i^x,*,^ > d t 6ipdxdt : 
Jt x Jo 

i-T2 poo 

■h = - / / F(u e )6d x ipdxdt 

-/ / < v x , t , F > 6 d x ip dxdt - 0(0) / (< /j yut ,F > -d y < n,U >\ y = yi )6dt, 

JT X Jo Jt x 

pTz poo PT2 

I / U(u e )9d xx <pdxdt - d x <p(eyi) / U(v e )\ y=yi 9 dt 



and 

„ = - e 
-> 0. 
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Since Ej v < and 



pT% pOO pT?2 

I I < u Xjt ,F > 9d x (pdxdt = / <v ,F > 9dt + 0(l)\\<p\\ L i, 
Jt x Jo Jt x 

we obtain an inequality of the form 

"T2 pT<2 

(< lh u t,F> -d y <fji,U >\ y=yi )0dt > <p(0) \ 

'Ti JTx 



rp rp 

0(0) [ 2 (<^y 1 ,t,F>-d y < f x,U>\ y=yi )9dt > <p{0) [ 2 <u ,F> 6dt + 0(l)\\<p\\ L i, 

(2.30) 

which proves the third inequality in (2.8) by chosing dp > such that H^Hl 1 —> but </3(0) > 0. 
This complete the proof of Theorem 2.1. □ 

Remark. Additional uniform estimates and regularity can be obtained from the identity in 
Step 2 of the proof of Theorem 2.1. Let (U, F) be a non-negative entropy pair that is uniformly 
convex on U. Use the identity (2.23)-(2.24) with 

= 1, T 1 = 0, T 2 = T, <p=\, yi = 0, y 2 = 00. 

We assume additonally here that, for a fixed state -u^ and for all t, 

u e (x, t) — > ttoo, u x (x, t) — > as x — > 00. 

The initial data uj should also decay rapidly at infinity. We obtain the following identity 

T poc pT 

e I U(v £ (y,T))dy-e / U(v e (y,0)) dy + / F( Uoo )dt 



pi poo pi 

/ U(v e (y,T))dy-e / U(v*(y, 0)) dy + / F( Uoo ) 
Jo Jo Jo 

pT pT poo 

- (F(u B ) -VU(u B )d y v e (0,.)dt + / V 2 U(v e ) ■ (dyv e ,d y v e )dydt = 
Jo Jo Jo 

Since the following two terms are uniformly bounded 

poo poo 

|e / U(v e (y,0))dy\ = |e / U(uj)dx \ < 0(1), 
JO Jo 



T 

W(?4)<V; e (0,.)cft| < 0(1), 



(Cf.(2.26) with 6 = 1), we deduce the uniform bounds 



e 



/ U(v e (T))dy+ / V 2 £/(t> £ ) • (dyv e ,d y v e ) dydt < 0(1). (2.31) 
Jo Jo Jo 
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For every Lipschitz continuous function g, it follows from (2.31) that the sequence d y g(y e ) 
is bounded in L 2 , so converges weakly to a limit which is nothing but d y < g >: 

d y g(v e ) -> d y < /j,g> weak in L 2 (R x R + ). (2.32) 

□ 

2.2. Finite Difference Schemes. We now extend the above analysis to several classes of 
finite difference schemes that are known to be consistent with the entropy inequality (1.3). 
Theorem 2.3 below deals with the entropy flux-splittings introduced by Chen-LeFloch [9], 
which also includes as a special case the Lax Friedrichs type schemes. We treat the Godunov 
scheme in Theorem 2.4. 

We are given two mesh parameters r and h with A = r/h kept constant and small enough 
in order to garantee the stability of the scheme. We define the approximate solutions u h (x,t) 
by the scheme 

u h (x, t + r) = u h (x, t) - \g(u h (x, t),u h (x + h, t)) + Xg(u h (x -h,t), u h (x, £)) (2.33) 

and the initial and boundary conditions: 

u h (x,t) = ut(x) for all t < r, 

, (2.34) 
u (x,t) = usit) for all x < h. 

By convention, the functions u h are right continuous. For the Lax-Friedrichs type schemes, 
the numerical flux g is given by 

SLax(^, w) = \{f{v) + f(w)) - Q(w - v), (2.35) 

where Q G (0, 1) is called the numerical coefficient of the scheme. (Symmetric positive definite 
matrices Q could also be dealt with.) For the flux-splitting schemes, g takes the form 

^piit(^,^) = /"H + / + (^), (2.36) 

where / = /~ + / + is a given entropy flux-splitting for the system (1.1). By definition [9], 
the matrix V/^ have real eignevalues and a basis of eigenvectors and there exists a pair of 
functions F± such that (U, F^) is an entropy pair for the system associated with flux- functions 
f ± . Observe that (2.35) is a special case of (2.36) as was pointed out by Chen-LeFloch. 
As in the analysis of Section 2.1, we assume a uniform L°° bound: 

\\u h \\ L ~(RxR + ) <0(1). (2.37) 

We rescale u h and define the function v h : R x R + — > U by 

v h (y,t) =u h (yh,t) y>0, t>0. 
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Let v and n be two Young measures associated with u h and v h , respectively. 

The entropy flux-splitting schemes satisfy discrete entropy inequalities of the form 

U(u h (x, t + r))- U(u h (x, t + r)) + X (g{u h {x, t),u h (x + h,t))- G(u h (x - h, t),u h (x, t)fj < 0, 

(2.38) 

where G is called the numerical entropy flux. With obvious notation, we have 

G hax (v, w) = ±(F(v) + F(w)) - Q(U(w) - U(v)) (2.35bis) 

and 

Gsput {v,w) = F~ (w) + F+ (v) . (2.36bis) 

Note that (2.38) hold for (2.36)-(2.36bis) provided u takes its value in a sufficiently small 
neighborhood of a given state in U. This is in constrast with the vanishing viscosity method 
where no such assumption was necessary. 

Theorem 2.1 admits the following extension to the flux-splitting schemes. We omit the proof 
which follows the lines of the one of Theorem 2.1. 

Theorem 2.3. Assume thatlA is a small neighborhood of a constant state inR N . The measure 
Hy : t is defined for ally > and almost every t, and is constant for y G [k, k + 1) for any integer 
k. For all convex entropy pairs (U, F), all y > 0, and in the sense of distributions in t G R+, 
one has 

F + (U B )+ < fil,t, F~ >> < Hy, t , F + > + < fiy +1 , t , F~ > 

> < fjLy +1 , t , F+> + < fi y+2jt , F~ > (2.39) 

> < vo, t ,F>, 

W),t = S UB (t) f° r a - e - t> 0, (2.40) 

and 

lim < fJ, y ,t, F + > + < fJ> y +i,t, F~ >>< VQf, F > . (2.41) 

□ 

Consider next the Godunov scheme corresponding to the flux g given by 

#Godunov(f , w) = f{R{v, w)), (2.42) 

where we denote by R(v,w) the value at x/t = 0+ of the solution to the Riemann problem 
with v and w as left and right initial data, respectively. The entropy flux is 

G G odunov(f , w) = F(R(v, w)), (2.42bis) 

Here it is more convenient to consider the values R(u h (x, t), u h (x + h, t)) and define a function 

w h 

w h (y, t) = R(u h (yh, t), u h (yh + h,t)) (2.43) 

for all y > 0. We denote by 7r a Young measure associated with w h and by v a Young measure 
for u . It is not difficult to extend Theorem 2.3 as follows: 
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Theorem 2.4. The measure n y j is defined for ally > 1/2 and almost every t, and is constant 
in y for y <E [k — 1/2, k + 1/2) for any integer k > 1. For all convex entropy pairs (U, F), all 
y > 1/2, and in the sense of distributions in t G R+, one has 

< K 1/2 ,t,F>> < 7iy,uF> 

> < TT y+ljt ,F> (2.44) 

> < v 0jt ,F>, 

and, at y = 1/2 and y = oo, n satisfies 



< 7T 1/2 ,t,F>= \imR(u B ,v h (l,t)), (2.45) 



and 



lim < Tfy ti F >>< u 0jt , F > . (2.46) 

□ 

We conclude this section by giving the main conditions satisfied by the discrete boundary 
layer, which will be studied in the rest of this paper. 

Assuming in the results of Theorem 2.3 that \i is a Dirac mass, say fj, = S v , the discrete 
boundary layer equation associated with the scheme (2.33) takes the form: 



9(v(y - l),v(y)) - g(v(y), v(y + 1)) = for all y > 1, 
v(y) = u B , ye [0,1), 



(2.47) 



while the discrete boundary entropy inequality is 

G(u B ,v 1 )>F(u ), (2.48) 

where v\ plays the role of a parameter. Formally, Theorem 2.4 leads to the same equations 
(2.47)-(2.48) with flux and entropy-fluxes given by (2.42). 

3. Sets of Admissible Boundary Values. 

Based on the results in Section 2, we introduce in this section several sets which can be 
used to formulate the boundary condition. For every method of approximation considered in 
Section 2, we introduce two different sets of admissible boundary values: 

(1) One is based on the entropy inequalities, £ entropy (uB) and yields a boundary condition 
of the form (1.5). This boundary condition is rigorously satisfied by the limiting func- 
tion generated by a sequence of approximate solution, as was proven in Section 2. For 
arbitrary systems having few or even just one entropy, the set S entropy (uB) may be too 
large to lead to a well-posed problem; 

(2) Another set, £ layer (ws), is based on the boundary layer equation, which was obtained 
formally after the analysis in Section 2. This set is more adapted to deal with general 
systems and lead to a well-posed problem when the boundary is non characteristic. 
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In this section, we study the local structure of those sets; under certain assumptions, we can 
prove that the sets £ layer (uB) are manifolds with dimension equal to the number of negative 
wave speeds of the system (1.1). This ensures that the initial-boundary value problem is well 
posed if, for instance, the data are constant states (boundary Riemann problem) as can be 
seen by applying the theory in [36]. We recall that (1.1) is assumed to be strictly hyperbolic 
throughout this section and we denote by Xj(u) the N real and distinct eigenvalues of the 
matrix V/(w) and by £j(u) and rj(u) corresponding basis of left and right eigenvectors. 

3.1 Vanishing Viscosity Method. For the sake of generality, we consider 

d t u e + d x f(u e ) = e d x {B{u e )d x u e ) , x > 0, t > 0. (3.1) 

Theorem 2.1 could be partially extended to this case. We assume that the viscosity matrix 
B(u) depends smoothly upon its argument u and is positive. We consider entropies U that 
are B-convex in the sense that V 2 U(u)B(u) > for all u under consideration. The boundary 
layer equation here takes the form 

dyf(v) = dy(B(v)dyV) (3.2) 

and the boundary entropy inequalities have the same form (2.16) but now U must be S-convex. 

Following Dubois-LeFloch [15], we introduce a set based on the boundary entropy inequal- 
ities. From now on, the time-dependence may be omitted. 

Definition 3.1. Given ub G U, the set of admissible boundary values based on the entropy 
inequalities associated with the vanishing viscosity method (3.1) is 

CsStyM = NeW; for all B-convex (U,F), F(u B )+VU(u B )(f(u )-f(u B )) > F(u )}. 

(3.3) 
□ 

It is obvious that this set may be quite large when the system (1.1) only admits few entropies. 
For most systems (N > 3), this set is too large to be used to formulate the boundary condition. 
In any case, it is difficult to get information on its local structure at u B . For general systems, 
the following observation is immediate. Fix a state u B G U and suppose that for some p one 
has 

\{u B ) < < X p+ i(u B ) (3.4) 

and the basis rj(u) is a family of eigenvectors for B(u). Then the set obtained by formally 
plugging the expansion 

f(uo) ~ f(u B ) + V/(ub)(m - u B ) + V 2 f{u B ) ■ (u - u B , u - u B ), 

o (3.5) 
F(u ) w F(u B ) + VF(u B )(u - u B ) + V F(u B ) • {u - u B , u - u B ) 

in the definition of ^'^(ub) contains ub+ the span of rj(us),j = 1, ---P and is contained 
in a cone with vertex u B - Indeed the inequality under consideration in (3.3) then becomes 

V 2 U(u B ) Vf(u B ) (u -u b ,uq- u B ) < 0. (3.6) 
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Since U is a convex entropy, the eigenvalues of V/(ws) satisfy (3.4), and V 2 U(ub) is a positive 
definite matrix, our claim follows. 

We also consider a second set of admissible boundary values, first introduced by Gisclon 
and Serre [20]. 

Definition 3.2. Given any ub G U, the set of admissible boundary values £ v igcosity ( u b)> based 
on the boundary layer equation associated with the vanishing viscosity method is the set of all 
Voo G U such that the problem 

B(v)d y v = f(v) - f(voo), 

v(0) = u B , (3.7) 
lim v(y) = Voo. 

admits a (smooth) solution v(y) G U for y > 0. □ 

To study the local structure of ^ v Tscosity(' UB )' we a PPly the following theorem concerning the 
existence of invariant manifolds. Cf. Hartman [23] for a proof. 



Theorem 3.3. Consider the differential equation 

££ + #(£,£o), H(y)eR N ,yeR, (3.8) 



dy 

where H : R N x R N — > R N is of class C 1 and for each £o 

#(0,£o) = ^(0,£o) = 0, (3.9) 

and E is a constant square matrix with d eigenvalues having negative real part, e eigenvalues 
having positive real part, and N — d — e eigenvalues having zero real part. For every (small 
enough) £o £ R N , let ^ y = £(y;£o) be the solution of (3.7) with the initial condition £(0;£o) = 
£o- Denote by T y the mapping £ — > £(j/;£o)- 

There exists a one-to-one mapping of class C 1 , S : £ — > S(£) = (w 1 , w 11 , w 111 ), having 
non-vanishing Jacobian and defined on a neighborhood of £ = G R onto a neighborhood 
of (w 1 ,w n ,w IH ) = (0,0,0) eR d x R N - d ~ e x R e , such that the mapping STyS' 1 takes the 
simple form 

STyS- 1 : w T y = e pI y W l + W\y- w 1 , w^ 1 , w IH ), 

^'"vw" 1 + W III {y-,w I Q ,w I Q I ,w I Q „ 

where P 1 , P 11 , and P 111 are constant real-valued matrices with all eigenvalues having moduli 
less than one so that the matrix exponentials e p , e p , and e p are well-defined, the absolute 



w y = er'^wi 1 + W II (y;w I ,w I I ,w I n ), (3.10) 
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value of any eigenvalue of e p is less than 1, and that for e p is greater than 1, and that for 
e p is exactly 1. Moreover the mapping W 1 , W n , and W IH are of class C 1 and their first 
order partial derivatives with respect to (w 1 , wig 7 , Wq 11 ) vanish at (0,0,0). Moreover one has 

W 1 = and W n = if w z = and w 1 ^ = 0, (3.11) 

and 

W n = and W 111 = if w^ 1 = and w^ 11 = 0. (3.12) 

□ 

The condition (3.11) means that the e-dimensional plane {wq = O.Wq 1 = 0} is a locally 
invariant manifold. If S(£o) belongs to this plane, then |£(y;£o)| ~^ 00 as V ~^ °°- The manifold 
{£ / Wq — 0, Wq 1 = 0} is called the unstable manifold of initial data for the equation (3.8). 

The condition (3.12) means that the <i-dimensional plane {wq 1 = 0,Wq H = 0} is a locally 
invariant manifold. If S(£q) belongs to this plane, then £(y;£o) —> as y — > oo. The manifold 

I w o T = 0' w q H = 0} i s called the stable manifold. 

Using Theorem 3.3 we prove the following result. 

Theorem 3.4. Let ub be given and assume that, for all u in a small neighborhood of ub, 

the basis rj{u) is a family of eigenvectors for B(u), 
the eigenvalues of B(u), saybj(u), are positive, ^ ^ 

and 

X P (u) < < X p+ i(u) (3.14) 

holds for some p. Then the set £viscosity( Ms ) contains the point ub and, locally nearby ub, 
contains a manifold with dimension p at least. When < \ p +i(u B ), ^v^cos\tyi UB ) ^ s a man tf°ld 
with dimension exactly p and its tangent space at the point ub is spanned by the eigenvectors 
rj(u B ), j = 1,2, ••• ,p. □ 

Assumption (3.13) holds for the examples considered later in this paper, but could easily 
be relaxed. A result similar to our Theorem 3.4 is also proved by Gisclon in [19], by another 
method. 

Proof of Theorem 3.4. The system in (3.6) can be written in the form 

— = Biy^VfMv + G(v, Voo ), 
ay 

v(0) = u B - Voo, (3.15) 
v(oo) = 0, 

where v(y) = v(y) — and the mapping G(v,Voo) satisfies G(0, Wqo) = 0, ^2(0,i>oo) = 0. In 
view of the assumption (3.13), the two matrices V/(v 00 ) and S(f 00 ) _1 V/(f 00 ) have the same 
eigenvectors, and so exactly the same number of positive, zero, and negative eigenvalues. Let 
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be the eigenvalues of -B(i>oo) 1 V/(v 00 ). Applying Theorem 3.3 with 

we see that there exists a one-to-one C 1 mapping 5, defined on a neighborhood of G R N , 
onto a neighborhood of (w 1 ,w H , w 111 ) = (0, 0, 0) G -R p xR N ~ p ~ 1 xR 1 , such that the manifold 

£ = {i;/u; JI (i;) = 0, w m (*5) = 0}, 

which is of dimension p, is stable. For any point ub — f oo taken in this manifold as an initial 
data for the differential equation in (3.15), the solution v(y) converges to as y — > oo, which 
is the third condition required in (3.15). 

If Woo belongs to this manifold, then (3.15) has a solution and hence solves the boundary 
layer problem. Furthermore the local structure of the set nearby ug can be described as follows. 

Suppose that < X p +i(ub)- The following estimate follows from (3.15): 

JV 

v(y) =Y j e Xi H j {v 00 ) ■ (u B - v^rjM + 0(v(y)) 2 . (3.16) 

For the right handside of (3.16) to go to zero, we must have 

9j(y oo ) = £ j (v oo )-(uB-v oo ) = 0, j =p+l r --N. (3.17) 

Keeping ub fixed, consider the map g : U — > R N ~ P with components gj given by (3.17). We 
have 

dg 

- — (u B ) = -(£ p +i(u B ), ■ ■ ■ ,£n(ub)), (3.18) 

(IV qo 

whose rank is N — p. By the implicit function theorem, (3.17) defines a manifold passing 
through ub and of dimension p. By construction its tangent space at wg coincides with the 
one for the stable manifold S. Therefore, in view of (3.18), the tangent space at u B for £ is 
spanned by the rj(u B ), j = 1, 2, • • • ,p. □ 

A general inclusion can be proven regarding the sets introduced in the previous sections. It 
has been first pointed out by Serre [44] (cf. also [19]) that: 

Proposition 3.5. The two family of sets introduced in Definitions 3.1 and 3.2 satisfy the 
inclusion 

^fcoWM C £SSty(«B) ( 3 - 19 ) 
for all ub EU. □ 
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Proof of Proposition 3.5. Let be a point in ^viscosity ('"-b) and denote by y — > v(y) the 
associated boundary layer function which satisfies v(0) = u B and v (oo) = v^. Consider the 
following function of the variable y > 0: 

n(y) = F( Voo ) - F(v(y)) - VU(v(y)) (f( Voo ) - f(v(y))). 

It is easy to see that 

^(y) = V 2 ^(y)) (/( V<x> ) - /( W (y)), fM - f(v(y)) 
> 

So the function O is non-decreasing, and since lim^^ O(y) = 0, we deduce that O(y) < for 
all y, in particular for y = 0, that is 

F(voo) - F(u B ) - VU(u B ){f( Voo ) - f(u B )) < 0. 

Thus Voo belongs to CX^)" D 

3.2 Finite Difference Schemes. We now turn to formulations of the boundary condition 
that are based on finite difference approximations. We use the notation in Section 2.2. We 
consider a scheme characterized by its mesh parameters r and h with A = r/h small enough, 
and by its numerical flux g(.,.) and its family of numerical entropy fluxes G(., .). It is tacitly 
assumed that the values u remain in a small neighborhood of a given state and attention is 
restricted to those entropies U such that the discrete entropy inequalities (2.38) are satisfied. 
In fact attention is mostly restricted to the Lax-Friedrichs type schemes and the Godunov 
scheme. 

Definition 3.6. Given u B G U, the set of admissible boundary values based on the entropy 
inequalities associated with difference scheme is 

£schemI( u B) = {^0 £ U\ there exists v\ s.t. for all convex (U,F), G(u B ,v\)) > F(u )}. 

(3.20) 
□ 

As for ^scheme ( Ub )i ^ ms se ^ ma y ^ e l ar S e to garantee that the boundary value problem 
is well posed. We also use the obvious notation ^L^ ropy (ws), ^pnttin^^s)' an< ^ ^Godunov ( u b)- 

For general systems and the diagonalizable splittings, i.e. those such that the vectors rj 
form a basis of eigenvectors for the matrices V/ ± , we have the following fact. Consider a 
Lax-Friedrichs type scheme or, more generally a diagonalizable, entropy flux-splitting scheme. 
Fix a state u B G U and suppose that (3.4) holds for some p. Then the set obtained by 
formally linearizing the inequalities in the definition of ^scheme ( Ub ) contains w^+the span of 
rj(u B ),j = 1, ...p and is contained in a cone with the vertex at u B .To see this we formally plug 
the second order expansion 

F ± (u ) w F ± {u B ) + VF ± (u B )(u - u B ) + V 2 F ± (u b )(u - u B ,u - u B ) (3.21) 
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and obtain the second order version of the inequalities in (3.20): 

VF(« B )(«o-«B)+V 2 F(tt B ) (u -u B ,u -u B ) < V F~ (ub)(vi-u b )+V 2 F~ (u B ) (vi-u B ,vi-u B )- 

Using the trivial entropies (i.e. choose for F the components of /), we get an (second order) 
expression for v±: 

V "(w B )(fi-w B )+V 2 f~ '(ub)(vi-u b ,v 1 -u b ) = ^ f{u B ){uo-u B )+V 2 f(u B )(uo-u B -: Uq-Ub), 
which can be used to rewrite the above inequality: 

V 2 U(u B )V 'f(u B )(u - u B ,u - u B ) < V 2 U(u B )'yf~(uB){vi-UB,vi-UB)- 
At the first order, v\ is given by 

V/~(ub)(vi -«b) = ^f(u B ){uo-u B ) 
so we arrive at the inequality 

-Vf + (uB) T Vf-(uB)- T V 2 U(uB)Vf(u B ){uo-UB,Uo-UB) < 0. 

The desired result follows immediatly since rj is a basis of eigenvectors for the matrices V/ + , 
V/~, and V/, and the function U is convex. 

The second family of sets is now defined. 

Definition 3.7. Given any ub G U, the set of admissible boundary values ^heme( u B)' based 
on the boundary layer equation associated with the difference scheme is the set of all G U 
such that the problem 

g(v(y),v(y + l)) = f(voo), 

v(y) = u B for ye [0,1), ( 32 2) 
lim v(y) = v^, 

y^oo 

admits a (piecewise constant) solution v(y) G U for y > 0. □ 

To study the local structure of Scheme ( u b), we apply the following theorem concerning the 
existence of discrete invariant manifolds. (Cf. Hartman [23] for a proof.) 

Theorem 3.8. Let T : R N — > R N , £ — > £i> be a mapping of the form 

^1 =rto + E(Zo), (3.23) 

where E(£q) is of class C 1 for small £o an d satisfy E(0) = and j^(0) = 0, and the matrix 
T is constant, non-singular, and has d > 0, N — d — e, e > eigenvalues of absolute value less 
than 1, equal to 1, and greater than 1, respectively. 
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and 



There exists a map S of a neighborhood of £o = onto a neighborhood of the origin in the 
space of (wq^Wq 1 jWq 11 ) G R d x n N - d ~ e x R e such that S is of class C 1 with non-vanishing 
Jacobian and STS~ X takes the simple form 

STS' 1 : w{ = A^i + W 1 (wl , wi 1 , wi H ) , 

w{ T = A'^i 1 + W 11 ^, wg, wl 11 ), (3.24) 
w[ n = A^wi 11 + W III (w I 0l w I I 1 w I 11 ), 

where P 1 , P 11 , and P 111 are d x d, (N — d — e) x (AT — d — e), and e x e square matrices 
with eigenvalues of absolute value less than 1, equal to 1, greater than 1, respectively, and the 
mapping W 1 , W n , and W 111 are of class C 1 and their first order partial derivatives with 
respect to (wq, Wq 1 : Wq 11 ) vanish at (0,0,0). Moreover one has 

W 1 = and W n = if = and w 1 ^ = 0, (3.25) 

W H = and W IH = if w^ 1 = and w^ 1 = 0. (3.26) 

□ 

The condition (3.25) means that the plane vq = 0, wq = of dimension d is locally invariant 
manifold and if R(£o) belongs to this manifold then T n £o — > as n — > oo. 

The condition (3.26) means that the plane uq = 0, wq = is a locally invariant manifold 
and if R(£o) belongs to this manifold, | T n £o |^ oo as n — > oo. 

Using this theorem we shall prove: 

Theorem 3.9. Consider a Lax-Friedrichs type scheme. Let ub be given and assume that 
(3.14) holds for some p. Then the set ^1^{ub) contains the point ub and, locally nearby 
ub, contains a manifold with dimension p. When < \ p +i(ub), ^l1^{ u b) is a manifold 
with dimension exactly p and its tangent space at the point ub is spanned by the eigenvectors 
rj(u B ), j = 1,2,--- ,p. □ 

Proof of Theorem 3.9. We search for all that solve the problem: 

H(v(y),v(y+l), Voo ) = 

v(0) = 0, (3.27) 

f (OO) = Voo 



with 



H(v(y),v(y + 1), Voo ) = v(y + 1) - v(y) - A (f( v (y)) + f(v(y + 1)) - 2f( Voo )). (3.28) 
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Using the notation H = H(v, w, t>oo), we compute 

dH A 

-£-(v,W,Voo) =Id+ 7^W(», 

dH A (3 - 29) 

— (v,w, Voo )=Id-— Wf(w). 

For A/ (2Q) small enough, the matrix dH/dw is invertible and its inverse is uniformly bounded 
w.r.t the variables v, w, and v^. By the global implicit function theorem (see J.T. Schwartz 
[46]) the system (3.27) can be solved for v(y+l). So there exists a smooth mapping K(v(y), v^) 
such that 

v(y+l) = K(v(y), Voo ) (3.30) 
and K^oo^Voo) = 0. Moreover one has 

0K (v(y), <U = (id - ±Vf(v(y + l)))" 1 (Id + ±Vf(v(y))). (3.31) 



dv v K J ' v 2Q Jy v " "' v 2Q 
The system (3.30) can be linearized around v^: 

v(v+1) = . (Id -^ Jv ^ (Id+ ^ MHv) 

Set w*(y + 1) = v(y + 1) — foo. The system can be written as 

v - (y + i) = -(/rf- ^ v (0oo)) -. (H+ ^| (lUK(!/) 



In other words 
where 

and 

We observe that 



v*(y + 1) = A(v 00 )v*(y) + K*{v*{y), Voo ), (3.32) 

A( Voo ) ee (ld-±Vf(v oc ))- 1 (ld+ T^^M) (3.33a) 
OK* 

K* and - — -- vanish at v*(y) = 0. (3.33b) 
ov*(y) 



the eigenvalues of the matrix A(voo) are ^ — -, (3.34) 



1+ ^K(Voo) 
1 - ^Ai(Voo)' 
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where A^i^) are the eigenvalues of V/(v 00 ). 

Namely (3.34) follows from the fact that the following two statements are equivalent : 

(1) a is an eigenvalue of A(voo); 

(2) there exists r^O such that A(v 00 )r = ar. 

Using the expression (3.33a) of A(voo) and simplifying the resulting equation, we get 

A(l + a) 

So a is an eigenvalue of A(i' 0O ) if and only if ~^g~jp^y is an eigenvalue of dfduivoo) with right 
eigenvector r; so 

20(0-1) w , 

^nr = Ai(Voo) (3 - 35) 

for some z with left eigenvector liiy^) and right eigenvector rj(voo). Solving (3.35) for a we 
get ith eigenvalue of A(voo) 

i j 

2Q' 



1 + ^A^Voo) 



d; = 



2qW{Voo) 



(3.36) 



Let T be a matrix which diagonalize V/(t>oo). Then the same matrix diagonalize Aivoo): 

TAT~ X = diag(ai, a 2 , • • -a n ). 



Set w*(y + 1) = 7V(y + 1), we get 
w*(y+l) = 



\ 



a 2 



w*(y) + G*(T-'w*(y),v 00 ) 



where G* and " ' 
Note that 

and 



\0 a n / 



g^Tj are zero at = 0- 



ai < a2 < • • • a p < 1 < a p +i < ■ • • < a n . 



a p+ i = 1 <^> A p+ i(voo) = 0. 



(3.37) 



Since all the hypothesis of Theorem 3.8 are satisfied, there exists a p-dimensional invariant 
manifold defined near such that, if the data Vq belongs to this manifold, then w*(y + 1) — > 
as j/ — > oo. In fact in terms of the original variable v(y + 1), we have the expansion 



TV 



v(y + l)- Voo = J2 aV j < £ j( u )' Vb-Voo> rj ( Voo ) + 0(| v(y + 1) - Voc |) 2 . (3.38) 
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In order for this to go to zero, as y — > we must have 

<tj(v oo ),u B -v oo >=0, j=p+l,---N. (3.39) 

This for fixed ub defines a map from R N R N ~ P and whose Jacobian at ub = t>oo is the 
matrix whose N—p rows are £j(voo). Since £j(voo) are linearly independent by implicit function 
theorem we deduce that (3.39) defines a p dimensional manifold passing through ub and if 
is in this manifold then there exist a solution to (3.29) whose local structure is given by (3.39). 
□ 

The following general inclusion can be proven: 

Proposition 3.10. The two family of sets introduced in Definitions 3.6 and 3.7 satisfy, for 
all ub € U, 

S^LeM C CheTe 7 ^)- (3-40) 

□ 

Proof of Proposition 3.10. We consider as before a difference scheme that satisfies discrete 
entropy inequalities. For every Voo in the set Scheme (' Us ' there exists a corresponding boundary 
layer profile v(y), solution of 

9(v(y),v(y+l)) = f(voo). 
The function v(y) is actually a stationnary solution to the scheme since 

v(y) - v(y) + X(g(v(y),v(y + 1)) - g(v(y - l),v(y))) = 0. 

Therefore for every convex entropy pair (U, F), it satifies the entropy inequality 

U(v(y)) - U(v(y)) + \(G(v(y), v(y + 1)) - G(v(y - 1), v(y))) < 0, 

which is nothing but 

G(v(y),v(y+l))-G(v(y-l),v(y) < 
Since lim^oo v(y) =v OQ , we get 

G(v(y),v(y+1)) > F( Voo ) 
and so with y = 0, since v(y) = ub for y e [0, 1), 

G(u B , Vl )) > F(u ) 

with v\ = v(l). That establishes that belongs to the set £schemJ( u B)- ^ 

Finally we treat the Godunov scheme. The sets ^Godunov(' UB ) an< ^ ^Godunov^- 6 ) are defined 
by Definitions 3.6 and 3.7. We now prove: 
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Theorem 3.11. Consider the Godunov scheme and let ub £W be given. We have 

tt-ovM = (3-41) 

This set can also be described as the set 

^Riemann^ = [R^yj) / w eU }, 

where R(ub,w) denotes the value at x/t = 0+ of the solution of the Riemann problem with 
data ub and w on the left and right, respectively. Moreover when (3.4) holds for some p, the 
set above contains the point ub and, locally nearby ub, is a manifold with dimension p and 
with tangent space at the point ub spanned by the eigenvectors Vj(u B ), j = 1,2, • • • ,p. □ 

Observe that the Godunov scheme does not produce any boundary layer, in the sense that 
the layer contains no interior point. 

Proof of Theorem 3.11. We recall that the set ^Godunov( UB ) is defined by the equation 

f(u B ) = f(R(v(y),v(y+l)), 

v{y) = u B for all y G [0, 1), ( 3 _ 42 ) 
lim v(y) = Woo, 

while the set ^Godunov ( u b) is defined by the inequalities 

F(R(u B , vi)) > F(u ) for all convex pair (U, F) (3.43) 
and for some v\ G U. So it is not hard to see from the definition that 

S^—(U B ) C ^dunov(^). 

On the other hand the inclusion 

also holds in view of Proposition 3.10. 
It remains to show that 

CoZLK) c £ Riemann (U B ). 

Consider a pair (uo, ^i) that solves (3.43). Then we need show that there exists w such that 

R(u B , w) = u Q . (3.44) 

Using the trivial entropies, we get 

f{R{u B ,v 1 )) = f(u ) 

which, combined with the inequality (3.43), shows that the pair of states (R(u B , vi), uo) is an 
entropy satisfying, stationary shock wave. On the other hand the Riemann problem with left 
and right initial data u B and R(ub,v±), respectively, contains only waves with non-positive 
speeds. Therefore the Riemann solution, with «b as a left state and uq as a right state, 
only contains waves with non-positive speeds. This function takes the value uq in the whole 
half-interval x/t > and thus R(u B , uo) = u o, which proves (3.44) with w = uq. □ 
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4. Selected Examples. 

In this section, we consider first the convex scalar conservation laws and establish that all 
the sets introduced in Section 2 are essentially the same. Some remarks are then given for the 
linear hyperbolic systems. Next we return to the scalar equation and treat a non-convex flux 
function, showing again that the sets are the same with the exception of the set based on the 
boundary layer equations. Finally we treat the elastodynamics system and isentropic Euler 
system. 

4.1. Scalar Conservation laws: Convex Fluxes. We consider a scalar conservation law 
with strictly convex flux, i.e. / (u) > and analyze the boundary layer equation. Let u* 
be the unique point such that /'(«*) = 0. To the state ub, when ub ^ w*, we associate the 
solution u* B 7^ ub of the equation f{u* B ) = f(us)- 

We state here a theorem which says that some of the sets introduced in Section 3 coincide 
in this case. We also recover the formulation of the boundary condition discovered by Bardos- 
Leroux-Nedelec [4] and Leroux [35]. 

Theorem 4.1. Consider a scalar conservation laws with convex flux. 

1 ) For any ub = R, the sets of admissible boundary values ^v£cos^ty(' u - B )' ^Godunov ( u b)> 
and ^Godunov (^s) > coincide with 



'Riemann 



(-oo, u* B ] U {u B } i/«s>'U*, 
( oo, it*] if u B < m*. 



and 

Viscosity ( U B) = Niemann («b) 

2) Givenus G U = [— M, M] for a fixed value of M > 0, we set ||/'||<x> = su P w e[-8M,8M] \f'( w )\ 
and consider a Lax-Friedrichs type scheme with coefficient A and Q satisfying ||/'||oo<VQ — ^> 
then 

4ax r («s) n [-M, M] = £ Rieman > B ) n [-M, M] \ {u* B } 
£lf™ py (u B ) n [-M, M] = £ Riemann (u B ) n [-M, M] 

□ 

4.2 Linear Hyperbolic Systems. 

It is not hard to prove that for a linear and strictly hyperbolic system, the sets defined in 
Section 3 are all equivalent when boundary is not characterestic. We only consider here the 
case of the discrete boundary layer based on the Lax-Friedrichs scheme. 

We also focus attention in this section to establish that the restriction (3.12) on the viscosity 
matrix is essential to our purpose here, as was observed in another context by Majda-Pego [41] 
in their study of traveling wave solutions to (2.1). The following example shows a situation 
where the viscosity matrix is a positive diagonal matrix, and does not satisfy (3.12), while the 
formulation may lead to a "wrong" boundary condition. 
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We consider the linear system 



5 5\ /5 



d t u+[ _ 3 3 Jd x u = el Q 1 )d xx u. (4.1) 



According to our earlier analysis, the boundary layer equation is 

1/5 0\ / - 5 5 N 



/1/5 0\ / -5 5\ 
=1 1 ) I -3 3) d v v (y^ 



1 * e * 1 1 

0yi/*>(y)=(_3 ^9 v v(y). 

Integrating this equation once and using v(+oo) = Voo, we get 

My)=(l 3 3 ) ( 4 - 2 ) 

/ - 5 5\ 

Now the eigenvalues of I I are Ai = —2 and A2 = 0. On the other hand, the eigen- 

/ - 1 1 \ 

values of I J are jiii = and /i 2 = 2. The solution of (4.2) with the initial condition 

v(0) = v B - i>oo is 

v(y) - =< h,v B -Voc > fi+ < £ 2 , v B - Vqo > f 2 e 2y , 

where 

? ' = (t-0- 4 = (7171)' fl = (!lv1)' f2 = (s2 



In order for v(y) — > as j/ — > 00, we must have < £i,vb — v<x >= and < £2, — 1> 00 >= 
which means that = vb- This requires that we prescribe u at the boundary.This is wrong 
boundary condition for the hyperbolic system 

/ -5 5\ 
dtU + ( _ 3 3 ) ^ w = 

because none of the characterestics are entering. 

Let us now consider the numerical boundary layer for a general linear and strictly hyperbolic 
system. Set f(u) = Au, where A is a constant matrix. The boundary layer equation becomes 

^Av(y+l) + ^Av(y)-^(v(y+l)-v(y)) = \Av 00 , (4.3) 
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v(0) = V B , v(oo) = Voo. 

For a given ub, we search for the set of states for which this problem has a solution. Set 
v°°{y) = v(y + 1) — Uqo. The first equation in (4.3) becomes 

(XA - I)v°°(y) = -(XA + I)v°°(y - 1). (4.4) 

Let £j and rj be the left- and right- eigenvectors for A associated with the eigenvalues Xj. Set 
C j (y) =< £j,v(y + 1) >. From (4.4) we get 

(l-WJCiy)* = (l + \\ j )C*(y-l) 



or 



with 

Integrating this, we get 

C j (y) =< e^VB-Voo > Q*^ J ) 

or 

+ 1) - = v°°(y) = ( i^ av ) < ^' WB ~ Uo °^' 

For u(y + 1) — > u^, we need < £j,vs — foo >= 0, j = p+ 1, • • • n because Ai < A2 < • • • A p < 
< Ap + i < • • • A n . This gives correct boundary condition when the eigenvalues are not zero; 
i.e. to prescribe 

<£j,u> for j = p + 1, • • • , N. 

4.3 Scalar Conservation Laws: Non-Convex Fluxes. 

We return to scalar conservation laws but now with non-convex fluxes. For definiteness we 
treat the case of the cubic flux given by 

f(u) = ±(u 3 -3u), (4.5) 
which has one minima and one maxima; indeed 

/(l) = -1, f'(l) = 0, /"(l) = 3, /(-l) = 1, f'(-l) = 0, /"(-I) = -3- 
For a given ub G R and the function / given by (4.5), we shall need the solution of the equation 

f(u) = f(u B ),u^u B . (4.6) 
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If u B < -2 or u B > 2, there is no solution for (4.6). If u B G (-2, -1) U (1, 2), then (4.6) has 
exactly two solutions. In this case we denote by u B and u s B the largest and smallest solutions 
of (4.6), respectively. If u B = —2,-1,1, or 2, then (4.6) has exactly one solution; namely 
1,2,-2, and —1, respectively. 

For the formulation of the results in this subsection, it will be convenient to introduce the 
following set, which is either the empty set or contains a single element: 

0, if u B G (-oo,-2) U [-1,1] U (2, oo) 

{1} iiu B = -2 

E(u B ) = { {u s B }, if-2<u B <-l (4.7) 

{u B }, Hl<u B <2. 

{-1} Hub = 2. 

Theorem 4.2. Consider the scalar conservation law with the non- convex flux (4.5). 

1) For any u B G hi — R, the set of admissible boundary values ^viscosityO^)' ^Godunov^- 6 )' 
and ^Godunov^- 6 ) coincide with 

{wb}, 
{-2,1}, 
[u s B ,l]u{u B }, 
(u B )= { [-1,1] 

[-l,u B ]u{u B }, 
{ub}, 

. I 2 '- 1 }' 



'Riemann 



«/«B < -2 
= —2 
«/- 2 <u B < -1 
if - 1 < u B < 1 
if 1 < u B < 2 
if ub > 2 
if u B = 2 



and 



S, 



Layer 
viscosity 



(u B ) = £ 



Riemann 



(u B ) - E(u B )- 



2) Given any state u B <EU = [-M, M] for a fixed value M > 2, we set || /' ||oo= sn P w e[-SM,8M] 
f'(w) | and consider a Lax-Friedrichs type scheme with coefficient A and Q satisfying || /' ||oo 



^ < 1. Then 



4ax r («s) n [-M, M] = £ Rieman > B ) n [-M, M] \ E{u B ), 
C op >b) n [-M, M] = S Riemann (u B ) n [-M, M]. 



□ 
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The proof of this is straightforward and is omitted. 
4.4 Nonlinear Elastodynamics. The system considered now arises in the modeling of 
elastic materials [10]: 

d t v — d x u = 0, 
d t u - d x a{v) = 0. (4 ' 8) 

It describes the evolution of a nonlinear material with deformation gradient v and velocity u. 
The stress function a is assumed to be smooth enough and satisfy the following conditions: 

<x»>0, va"(v)>0. (4.9) 

Let us discuss the vanishing viscosity approximation for the viscosity matrix B(u) = I. The 
boundary layer problem to be studied here is 

— d y u = dyV, 
-d y a(v) = d 2 y u, (4iQ) 

v(0) = V B , f (OO) = Woo, 

u(0) = u B , u(oo) = Woo 

We need determine the set of (i>oo> w oo) f° r which (4.10) has a solution. Integrating once the 
ODE'S and using the boundary condition at infinity, we get 

d y v = u 00 -u, Uy = a(voo) - a(v). (4.11) 

Cross multiplying the equations and integrating, we get 

(it - Woo) 2 



I {<r{8)-<T{v 00 ))d8, 



so 

(u- Uoo ) = ±Q^ 2(<r(s)-<T(v 00 ))ds y j . (4.12) 

Note that f y (cr(s) — a^^^ds > because of the condition cr'(v) > 0. From (4.12) it follows 
that 

v(y) = v 00 ^ u(y)=u 00 . 

Since we are interested in a solution connecting (vb,ub) at y = to (foo 5 ^oo) at y = oo, we 
get from (4.11) that either 

vb < v(y) < and u B < u(y) < 

or (4.13) 

vb > v(y) > and u B > u(y) > u^. 
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This determines the sign in (4.12): 

/ rv \ 1/2 

Uoo + ( / 2 (a(s) - a(v)) ds J ift> > d c 



w = < 



u<x — I 2 (cr(s) - cr(v)) ds 

\J Woo 



1/2 



iff < w c 



Since we need (vb,ub) to be on this curve, we obtain that the set of (foo,Woo) so that (4.10) 
has a solution lies on the curve 



v B \ 1/2 

iii; - | / 2 (cr(s) - cr(w)) (is if u oo < u B 



u B + { 2 (cr(s) - cr(w)) ds J if ^ > w B , 



(4.14) 



where (vb,ub) is fixed. 



Let us now turn to the Lax Friedrichs scheme. For the system (4.8), the discrete boundary 
layer equation is 



H(v(y),v(y+ l),tioo,Voo) = 



A(w(w + 1) + w(w)) + u(w + 1) - - 2 Awoo) 
A(<r(v(y+1)) + (r(w(w)) + w(w+l) - u(y) -2Xa(v 00 ) 



(V,U)(0) = (Wfi, W ), (v,u)(oo) = (UocUoo). 

Here the eigenvalues of (4.8) are 

A 2 (Voo, Woo) = -Al(Voo,Woo) = Cr'(^oo) 1/2 

and, with the notations of Section 3, 

. . l-Xa'M 1 ' 2 . . 1 + Ao>oo) 1/2 

ail^oo, Woo) = — -— 77 TT77, 02(^00,^00) = 



= 0, 

(4.15) 



1 + Aa'^oo) 1 / 2 ' 



l-Aa'^oo) 1 / 2 ' 



Thus < 01(^00,^00) < 1, 02(^00, Woo) > 1- By the analysis of Section 3, it follows that the 
set of ( Wqo, Woo) near (v b ,ub) for which (4.15) has a solution lie on a curve passing through 
{v B ,u B )- 

4.5 Eulerian Isentropic Gas Dynamics. We now consider the isentropic approximation 
to the compressible Euler system. The system is composed of the two conservation laws for 
the mass and the momentum of a gas [10]: 



d t p + d x (pu) = 0, 

d t (pu) + d x (pu 2 +p(p)) = 0. 



(4.16) 
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The main unknowns are the specific density p and the velocity u. The pressure is a function 
of the density and, for simplicity, we shall restrict to a polytropic perfect gas: 

p(p)=p\ 7e(l,oo). (4.17) 

We consider the boundary layer equation generated by the vanishing viscosity method with 
B{u) = 1: 

d y (pu) = d 2 y p 

d y (pu 2 + p(p)) = d 2 u (4.18) 
p(0) = Pb, u(0) = u b , p(oo) = poo, u(oo) = u 00 . 

Integrating the ODE'S and using the boundary condition at infinity, we get 

d y p = pu - Poo Uoo 

d y u = pu 2 + p(p) - pooU 2 ^ - pipoo) (4.19) 

p{0) = pB, U(0)=U B , P(OO) = poo, u(oo) = Uoo. 

The eigenvalues of the matrix obtained by linearizing the R.H.S. of (4.19) around (pocn^oo) 
are 

Al(poo, Uoo) = Uoo - c(poo), A 2 (/9oo, «oo) = Uoo + c(poo) (4.20) 

where c 2 (poo) = p'(p)- We have to distinguish between five different cases. We define the 
following regions in (p, w)-plane: 

Q/ = {(p, u) : u - c(p) < 0, u + c(p) < 0} 
= {(p,u) :u - c(p) < 0, u + c(p) = 0} 
Mm = {(p,u):u- c(p) <0,u + c(p) > 0} (4.21) 
^iv = {(p,u) : u - c(p) = 0, u + c(p) > 0} 
Oy = {(p,«):w- c(p) > 0, u + c(p) > 0} 

Thus in O/ both eigenvalues are negative, whereas in O77 one has Ai < 0, A 2 = 0. In fi///, 
one has Ai < 0, A 2 > 0, wheras in Q/\/, one has Ai = 0, A 2 > and in 1V 5 Ai > and A 2 > 0. 
Following the analysis that we did for the proof of Theorem 3.2, it is not hard to get the 
following local result. 

Case 1 : (pb,ub) G O/. In this case the set of (p C xj,w 00 ) close to (pb,ub) for which (4.19) 
has a solution is an open neighborhood of 

Case 2 : (pb,ub) £ ^n- In this case the set of (poo.Uoo) close to (pb,ub) for which (4.19) 
has a solution is a union of a two-dimensional region U in 0/ and a curve in Qui through 
(Pb,ub) intersecting U. 
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Case 3 : (pb,ub) G ^in- In this case the set of states (poo>^oo) close to (pb,ub) for which 
(4.19) has a solution is a curve through (p#, mb) 

Case 4 '■ (pbiUb) G H/y. In this case the set of states (pco^oo) near (pb,ub) for which 
(4.19) has a solution lies in a curve in Vim through This does not extend to Vl V - 

Case 5 : (pb,ub) G Oy. There cannot be any point (poojUoo) in Vly for which (4.19) has a 
solution. 

Next we consider the Lax-Friedrichs scheme. The discrete boundary layer problem to be 
solved is 

X(p(y + l)v(y + 1) + p{y)v{y)) - 2\ PooUoo - (p(y + 1) - p{y)) = 

X(p(y + l)u(y + l) 2 + p(y)u(y) 2 ) - 2Xp 00 u 2 00 - (p(y + l)v(y + 1) - p{y)v{y)) + X(p(p(y + l))+p(p(y)) 

(4.22) 

(Pb,u b ) given and (p,tt)(oo) = (poo,Woo)- 

Given (pb,ub) we determine (poo, - ^) close to (pb,ub) for which (4.22) has a solution. 
Following the analysis of the proof of Theorem 3.4, we get the eigenvalues of the linearized 
matrix at (poo^oo) are 



1 + AAi(poo, Woo) / \ 1 + AA 2 (Poo,«oo, 

ai = a^poojMooJ = — — -,a 2 = a 2 {p 00 ,u 00 ) = 



l-AAi(p ) " ' ' 1 - AA 2 (Poo, «oo) 

where Ai and A 2 are given by (4.20). If (poo,Woo) G fij,ai < l,a 2 < 1, if (poo^oo) £ 
0/7, ai < l,a 2 = 1, if (poc^oo) G 0///,ai < l,a 2 > 1, if (poo,Woo) G 0/y,ai = l,a 2 > 1 
and if (p, «oo) e Oy,ai > l,a 2 > 1. It follows from the proof of Theorem (3.4), that if 
(PBiUb) G Oj, then the set of states (poc^oo) near (pb,v>b) for which (4.22) has a solution 
connecting (ps,ws) to (pco^oo) is a neighborhood of (ps,ws). If (pb^ub) G O/j this set is 
a union of an open set £7 in 0/ and a curve in Vim through (pb,wb) which interset £7. If 
(PBiUb) £ O/jj this set of (pooj^oo) near (ps,ws) consists of a curve through (pbjUb) and 
if {pbiUb) G 0/y this set consists of a curve in through (ps,us). If (pb,ub) G Oy no 
point (pooi^oo) G Oy can be connected by a solution of (4.22) from (ps, its). 

4.6 Lagrangian Isentropic Gas Dynamics. Finally, we consider the system of gas 
dynamics in Lagrangian coordinates 

dm - c^m = 0, 

d t u + d x = 0, (4 ' 23) 
where u is the velocity and v > is the specific density. The eigenvalues of (4.23) are 

Ai = -- < 0, A 2 = - > 0; (4.24) 
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hence the boundary x = is not characteristic. 

The purpose of this section is to provide an explicit formula for the boundary layer set 
associated with the Lax-Friedrichs scheme. The boundary layer equation takes the form 

X(u(y + 1) + u{y)) -2\ Uoo + v(y + l)- v{y) = 

4t^tt + 4t)- 2 --^ +1 » + "(^° <426) 

\v(y+l) v(y)J Voo 

with 

(v(0), u(0)) = (v B , «b), (u,u)(oo) = (i; 00 ,u 00 ). (4.26) 

We restrict attention to vb > 8 > for fixed S, and we determine the set of (vqo, Moo) for which 
(4.25) has a solution. We set 

w(y) = ^ (4.27) 



so that (4.25) becomes 



1 2 

- u(y + 1) + u(y) 



w(y+l) w(y) ' Woo 

w{y + 1) - w(y) + u(y + 1) + u{y) = 2 u^. 
Adding the two equalities, we get 

w(y + l)+ 1 + - w(y) + 2u(y) = — +2u 

Setting 

N(y) = -2u(y) + 2 Uoo - — — + — + w(y), 

w{y) Woo 

we obtain a quadratic equation for w(y + 1): 

«; 2 (y + l)-Ar( 2 /)«;( 2 / + l) + l = 0. (4.28) 

Therefore 

W (y + l) = i(iV( 2/ )±(iV( ?/ ) 2 -4) 1/2 ) 
from which we get an expression for u(y + 1) as well: 

«(!/ + 1) = ^ J%) ± ^ (N(y) 2 - 4) 1/2 - (4.29) 
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Observe that iV(oo) = + l/Vx,, where = v^/X and N{oo) 2 — 4 = (w^ — l/w^) 2 . The 
product of the two roots of (4.28) is equal to one. Stability requires > 1 so we choose the 
larger root in (4.29). We have finally from (4.29) and (4.25). 

v(y+l) = ±N(y) + ±(N(y)*-4) 1 '* 

(4 30) 

u(y + 1) = 2 Uoo - u(y) + V -f - - \(N(yf - 4)^. 

The Jacobian of the R.H.S. of (4.30) at (uoo,i>oo) is easily seen to be 

A(Vooi Woo) = 



w 2 — 1 w 2 — 1 

oo oo 

-2 «& + 1 



whose eigenvalues are 
In terms of v^o, we have 



\A«-1) w^-lj 



Woo - 1 Woo + 1 

ai = — , a 2 = -. 

Woo + 1 Woo ~ 1 



1 — 1 + — 

ai ~ i + X' a2 = i _ X 



If the data for the Lax-Friedrichs scheme are chosen such that the v component is bounded 
away from zero, then so is the approximate solution. Hence we can restrict attnetion to Voo > S' 
for some 5' > 0. For A small enough, we have 

< ai < 1 and a 2 > 1, 

and Theorem 3.10 applies. We deduce that the set of all states (vqo, Uoo) near (vb,ub) for 
which (4.25)-(4.26) has a solution is a curve passing through (vb 7 ub)- 

5. Concluding Remarks. Given a family of sets such as those introduced in this paper, 
we can formulate the boundary condition for the hyperbolic problem. When the solutions u 
under consideration are functions of bounded variation, the traces exist in a strong sense and 
one can require that 

u(0+,t) E S(u B (t)), t>0, (5.1) 

holds for all, except countably many, t. This type of regularity has been recently rigorously 
established by Amadori in her thesis [1] , using the front tracking scheme and the sets £ Godunov 
(= ^godunov = ^Godunov)- ^ would be interesting to extend [1] to the other sets we introduced 
here. 

For general L°° solutions constructed by the vanishing viscosity method, the boundary 
condition 

supp v 0t t C CSty ( 5 - 2 ) 
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has been rigorously derived in Theorem 2.1. When the method of compensated compactness 
applies [12], an existence theorem for the boundary-value problem (1.1)-(1.3), (5.2) follows 
immediatly from Theorem 2.1. Such a result is satisfactory only when the condition (5.2) 
yields, for simple enough initial and boundary data at least, a well-posed problem. This is the 
case for the scalar equations and the linear systems, and, likely, for any system in the so-called 
Temple's class (having coinciding shock and rarefaction curves). 

In a recent preprint by Grenier and Gues, a scaling of the type 1/y/e is used for linear 
systems of equations to obtain a more precise description of the boundary features. However, 
as far as the formulation of a well-posed, limiting boundary-value problem for the hyperbolic 
equations is sought, the scaling 1/e we used in (2.5) happens to be sufficiently discriminating. 
The formulation based on the boundary layer equation may not be appropriate as it is when 
the boundary is characteristic. On the other hand, the formulation based on entropy inequality 
capture rigorously some of the features in the solution near the boundary, but is more difficult 
to work with analytically. Further study of the connection between the two sets for systems is 
in progress. 
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